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Abstract 

Simulation methods for the decays B — > tt + tt~V, where V is a 1 vector- 
meson, are presented in detail. Emphasis is put on the use of the helicity formalism 
and the use of effective Lagrangians. We show the importance of p° — to mixing in 
enhancing the direct CP violation (DCPV) when the pion-pion invariant mass is 
near the mass of the oj. 



1 Introduction 



In the framework of the LHCb experiment devoted to the search for CP violation 
and rare B decays, special care is given to the B decays into two vector mesons, 
B — > V1V2, Vi = 1 . Physical motivations for studying these processes are numerous: 

(i) Weak interaction governing the B decays, the vector-mesons are polarized and 
their final states have specific angular distributions; which allows one to cross-check the 
Standard Model (SM) predictions and to perform tests of models beyond the SM. 

(ii) In the special case of two neutral vector mesons with V° = V°, the orbital angular 
momentum, £, the total spin 5* and the CP eigenvalues are related by the following 
relations: 

e= s = 0,1,2 ^cp = {-if , 

which implies a mixing of different CP eigenstates, leading to CP non-conservation 
process. According to Dunietz et al [T], tests of CP violation in a model independent 
way can be performed and severe constraints on models beyond the SM can be set. 

2 Helicity formalism and its applications 

Because the B meson has spin 0, the final two vector mesons, V\ and V 2 , have the same 
helicity Ai = A2 = —1,0,-1-1, and their angular distribution is isotropic in the B rest 
frame. Let H w be the weak Hamiltonian which governs the B decays. Any transition 
amplitude between the initial and final states will have the following form: 

H X =(V 1 (\)V 2 (\)\H V> \B) , (1) 

where the common helicity is A = —1, 0, +1. Then, each vector meson V, will decay into 
two pseudo-scalar mesons, a^, 6j, where and hi can be either a pion or a kaon, and the 
angular distributions of a, and hi depend on the Vi polarization. 

The helicity frame of a vector- meson, Vi, is defined in the B rest frame such that the 
direction of the Z-axis is given by its momentum p^. Schematically, the whole process gets 
the form: 

B — ► V x + V 2 — ► (ai + 61) + (a 2 + h 2 ) . 

The corresponding decay amplitude, M\{B — > Y^=ii a i + i s factorized according to 
the relation, 

2 2 
M A B - X>i + h d) = H ^ B - v i + v *) x HMVi - <k + bi) , (2) 

8=1 8=1 

where the amplitudes Ai(Vi — > Oj + 6«) are related to the decay of the resonances Vi. The 
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Ai(Vi — > a,i + hi) are given by the following expressions: 

1 

A x (V 1 ^a x + b x ) = Yl c i^WOA,0), 

mi=— 1 
1 

A 2 (\/ 2 ^a 2 + 6 2 )= ^ c 2 i^ m2 (0,0 2 ,-0) . (3) 

r«2=— 1 

These equalities are an illustration of the Wigner-Eckart theorem. In Eq. (|5|L the Ci and 
c 2 coefficients represent, respectively, the dynamical decay parameters of the Vi and V2 
resonances. The term m . (0j, 0j, — <^) is the Wigner rotation matrix element for a spin-1 
particle and we let A(a,) and A(6j) be the respective helicities of the final particles Oj and 
6j in the rest frame. 61 is the polar angle of ai in the V\ helicity frame. The decay 
plane of V\ is identified with the (X-Z) plane and consequently the azimuthal angle <p\ is 
set to 0. Similarly, 62 and <p are respectively the polar and azimuthal angles of particle 
a 2 in the V2 helicity frame. Finally, the coefficients rrii are defined as: = A(eij) — A(&j) 
Our convention for the D\ m _(a , (3 , j) matrix element is given in Rose's book [2], namely: 

D\, mi {a,l3,l) = expH(Aa + m iT )] < mi (/3) . (4) 

The most general form of the decay amplitude Aii^B — > Y^=i( a i + M) ^ s a ^ n ear 
superposition of the previous amplitudes M\(B — > Ylt=i( a i + denoted by: 

2 2 
M(£?^^(a J + 6 4 )) =^M A (fi^^(a, + 6,)) . (5) 

4=1 A i=l 

The decay width, r(£? — > V1V2), can be computed by taking the square of the modulus, 



which involves the three kinematic parameters, 61,62 and 



This leads to the following general expression: 

1 2 2 

d 3 T(B^V 1 V 2 )<x\Y,M x (B^J2^ + b i)) =J2KyFx,y(0i)Gx,y(6 2 A) , (6) 



i=l A, A' 



which involves three density-matrices, h^x, F\ t \r(6i) an d ^ 



A, A' 



* The factor h\ t \i = H\H^, is an element of the density-matrix related to the B decay. 

* F\ t \'{6i) represents the density-matrix of the decay V\ — > ai + b\. 

* G*a,a'(^2, 0) represents the density-matrix of the decay V2 — > a 2 + b%. 

The analytic expression in Eq. © exhibits a very general form: it depends on neither 
the specific nature of the intermediate resonances nor their decay modes (except for the 
spin of the final particles). 

The previous calculations are illustrated by the reaction B° — > K*°p° where K*° — > 
K + 7i~ and p° — > 7r + 7r~. In this channel, since all the final particles have spin zero, the 
coefficients mi and m 2 , defined previously, are equal to zero. The three- fold differential 
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width has the following form: 
d 3 T(B -> V X V£ 



d(cos 8i)d(cos # 2 )c/0 
+ 



oc (h ++ + h )sin 2 6 l isin 2 6 l 2/4 + /i 00 cos 2 6 l 1 cos 2 ^ 2 



j<^fe(/i + o)cos — J%i(/i + o)sin0 + ^(/io-)cos0 — J%i(/i -)sin 0jsin 26 l isin 26> 2 /4 

+ j^O + _)cos 20 - J%i(/i + _)sin 20 jsin 2 # lS in 2 # 2 /2 • (7) 

It is worth noticing that the expression in Eq. (jJJ) is completely symmetric in 9\ and # 2 
and consequently, the angular distribution of a x in the V\ frame is identical to that of a 2 
in the frame. From Eq. (JJJ) the normalized probability distribution functions (pdf) of 
01, # 2 and can be derived and one finds: 

/(cos0i )2 ) = (3/i 00 - l)cos 2 6>i >2 + (1 - /i o) , 

#(0) = 1 + 2 ^(/i+_)cos 20-2 J%i(/i+_)sin 20 . (8) 



3 Final state interactions and p — oj mixing 

Hadrons produced from B decays are scattered again by their mutual strong interactions, 
which could modify completely their final wave-function. Computations of the branching 
ratios B — > Hadrons must take account of the final state interactions (FSI) |3| which are 
generally divided into two regimes: perturbative and non-perturbative. 

An important question arises: how to deal with the FSI in a simple and practical way 
in order to perform realistic and rigorous simulations? 

The method which has been followed for the simulations is largely developed in the Ref. |H] 
and is based on the hypothesis of Naive Factorization, which can be summarized as 
follows: 

• In the Feynman diagrams describing the B decays into hadrons like tree or penguin 
diagrams, the soft gluons exchanged among the quark lines are neglected. 

• Using the Effective Hamiltonian approach and applying the Operator product 
Expansion method (OPE), perturbative calculations are performed to the Wilson 
Coefficients (W.C.), Cj, at the Next to Leading Order (NLO) for an energy scale 
> m B . 

• Non-perturbative effects representing physical processes at an energy < m# are 
introduced through different form-factors. 

• The color number, N c , is no longer fixed and equal to 3. It is modified according to 
the following relation: 

1 1 

(iVF) " 3 + * ' 

where operator(s) £ describe(s) the non-factorizable effects. 
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Another important effect which appears in the channels B — > -k + -k^V is the p° — uo 
mixing which is an unavoidable quantum process. Indeed, the tree amplitude, A T , and 
the penguin, A p , are modified according to the following relations: 



(K*7r-n+\H T \B} = -^fl pw t UJ + ^t p} 

(K*n-n+\H p \B) = -^-ft^ + % p . (9) 



SpS^j Sp 



Here ty (V = p or uo) is the tree amplitude and py the penguin amplitude for producing 
a vector meson, V, g p is the coupling for p° — > 7r + 7r~ , 11^ is the effective p° — uo mixing 
amplitude, and sy is the inverse propagator of the vector meson V , sy = s—m v +imyTy 
with y^s being the invariant mass of the pair. 

The ratio A P /A T , which is a complex number, gets the final expression: 

re iS e i(t> = + s ^Pp 

where 5 is the total strong phase arising both from the p° — u resonance mixing and the 
penguin diagram quark loop and </> is the weak angle resulting from the CKM matrix 
elements. 



4 Explicit calculations, simulations and main results 

Computations of the matrix elements are based on the effective Hamiltonian given by: 

We// = ^ E VckmC^O^) , (11) 

* i 

where Gp is the Fermi constant, Vckm is the CKM matrix element, Ci{p) are the Wilson 
coefficients (W.C), Oi(p) are the operators related to the tree, penguin-QCD and penguin- 
EW diagrams and p represents the renormalization scale, which is taken equal to tub- 

The W.C are calculated perturbatively at NLO by renormalization group techniques 
[3], while the non-perturbative parts, related to the operators Oi and form factors. The 
latter are explicitly calculated in the framework of the pioneering BSW models [Sj. The 
free parameters which remain are: (i) the ratio q 2 /mt ) 2 , where q 2 is the invariant mass 
squared of the gluon appearing in the penguin diagram and (ii) the effective number of 
colors, N c eff . 

Combining both the Wilson coefficients and the BSW formalism by including the p°—uj 
mixing in the B meson rest-frame, where Pb = (ms,0), the helicity amplitude is given 
by the final expression: 

H X (B - p°(uj)V 2 ) = iB p x {V ub V: s c p tl - V tb V*cf pi ) + iC{{V ub V: s c p t2 - V tb V* s c^)+ 

iB" x {V ub V: s cl - V tb V*c» pi ) + iC»(V ub V: s cl - V tb V*c^ 



(s p - ml) + imj 



(12) 
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where the terms and C x * are combinations of different form factors. Their explicit 
expressions, corresponding to the helicity values (A = —1, 0, +1), are given in Ref. j^j. 

From the above expression, we can deduce the dynamical density-matrix elements, 
h\ x'-, which are given by: 

h x , x , = H X (B - p°(cu)V 2 )H* x ,(B - p°(cu)V 2 ) . 

Because of the hermiticity of the DM, only six elements need to be calculated. 

Main results 

* hi j elements depend essentially on the masses of the resonances; each resonance mass 
being generated according to a relativistic Breit- Wigner distribution: 

da_ =(J Y R M R 
dM* N (M*-Mlf + (T R M R f ' 

where 

• Cn is a normalization constant. 

• Mr and Tr are respectively the mass and the width of the vector meson. 

For computational reasons, the analytical treatment of p° — u mixing is simplified in 
the Monte-Carlo simulations [7j. 
-k The main conclusions are: 

1) The spectrum of hij is too wide because of the resonance widths, especially the p° 
width, T p = 150 MeV/c 2 '. 

2) The longitudinal polarization, h 00 = \H \ 2 , is largely dominant. 

In the case of B° — > p°(ui)K* , the mean value of h Q0 is ~ 87% while for B + — > p°(o;)p + , 
its mean value is ~ 90% . These results have been confirmed recently by both BaBar and 
Belle collaborations jH]. 

3) The matrix element h = |-^-i| 2 is very tiny, < 0.5%. 

4) The non-diagonal matrix elements hij are mainly characterized by: 

• The smallness of both their real and imaginary parts. 

• Jmj&fc w 0.001 -> 0.1. 

• In the special case of B + — > p°(u)p + , J%i(h;j) w 0.0. 

We arrive at the conclusion that there is a kind of universal behavior of the Density- 
Matrix Elements, whatever the decay B — > 7i + 7i~V is (V = i^*°,ii'* ± , p ± ). 
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Consequences for the angular distributions 

In the helicity frame of each vector- meson, Vt, the angular distributions given above (see 
Eq. (8)) become simplified: 

• According to the analytic expression for g(<p) and because of the small value of 
(h+-), the azimuthal angle distribution is rather fiat. 

• From the expression of /(cos^i^) and because of the dominant longitudinal part 
h 00 , the polar angle distribution is ~ cos 2 ^. 

Branching ratios and asymmetries 

* The energy E$ and the momentum pi of each vector meson vary significantly according 
to the generated event. So, the branching ratio of each channel must be computed by 
Monte-Carlo methods from the fundamental relation: 



dT(B -> V X V 2 



8vr 2 M 



\M(B -> V X V 2 



d 3 pi d 3 p 2 
2E\ ~2E 2 ~ 



5\P-p 1 -p 2 ) 



and 



■k For a fixed value of q 2 /rrib 2 , the BRs depend strongly on the Form Factor model. They 
could vary up to a factor 2. 

* The relative difference between two conjugate branching ratios, BriB — > /) and 
Br(B — > f), is almost independent of the form-factor models. 

* An interesting effect is found in the variation of the differential asymmetry with 
respect to the tttt invariant mass which is defined by: 

T m (B^f)-T m (B^f) 
aCPim) ~ T m (B^f)+T m (B^f) ■ 

acp{m) is amplified in the vicinity of the u resonance mass, a mass interval of 20 MeV/c 2 
around M w = 782MeV/c 2 . 

This differential asymmetry is ~ 15% in the case of B° — > K*°p°(uj) while it reaches 
90% in the channel B^ p ± p (u). It is almost independent of the form-factor model 
and the only explanation of this surprising effect is the mixing of the two resonances 
p° and oj . These results have already been predicted analytically in the channel B — > VP 
(B =>■ p°(u;)7r, B p°(u))K) by A.W. Thomas and his collaborators [Oj. 
Some physical consequences can be deduced: 

• A new method to detect and to measure the direct CP Violation both in B° and B ± 
decays can be exploited. 
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According to the analytical expression: 



dir A 2 — A 2 — 2 r sin 8 sini 

1 cp ~ 



A 2 + A 2 1 + r 2 + 2 r cos 5 cos 

an opportunity is offered for measuring sin $ , where $ is a weak angle resulting from the 
CKM matrix elements: 

★ $ = Arg[V tb V*/V ub V* s ] = 7 in the case of B -> p°(u)K*, 

★ $ = Ar#[l^yKbK*d = /5 + 7 = n - a in the case of 5 -> p°{uj)p ± . 

★ Another interesting result deduced from the above calculations is the ratio 
A P /A T [TU]. It depends on the free parameter q 2 jm 2 and almost independent (except in 
the vicinity of the uo resonance mass) of the tttt invariant mass: 

q 2 /m 2 = 0.3 < r > = 0.31 ± 0.03 , 
q 2 /m 2 b = 0.5 < r > = 0.27 ± 0.03 , 

while the standard estimation of the ratio A p / A T is ~ 30% (Buras et at). 



5 Comparison with recent experimental results 

The Belle and BaBar collaborations recently published their first results concerning the 
charmless B decays into vector mesons, B — ► VV |B|. 

Belle Collaboration 



Channel 


Br(xl0~ 6 ) 


/l = l^ol 2 


pV 


31.7 ±7. l(stat)t 3 6 8 7 (syst) 


0.948 ± 0.106(stat) ± 0.021(syst) 


Our results 


11.0 -> 20.0 


90% 



Babar Collaboration 



Channel 


Br(xl0~ 6 ) 


/x = l^ol 2 


A C p 


p°K*+ 


10.6i|^±2.4 


0.96i{j;$f±0.04 


0.20±g;|gd=0.04 


Our results 


2.3 -> 5.8 


87% 


-6.4% -> -22% 


pV 


22.5±|;I±5.8 


0.97l{j;{$±0.04 


-0.19±0.23±0.03 


Our results 


11.0 -> 20.0 


90% 


-8.5% -> -10% 
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Because there is as yet insufficient data to allow one to bin data in the region of the 
uj resonance, one can only look at the global asymmetry Aqp measured by BaBar. This 
is compatible with zero and the differential asymmetry with regard to rnv has not been 
investigated. We look forward with great anticipation to the time when the invariant 
mass distribution can be investigated. 

6 Conclusion and perspectives 

Monte-Carlo methods based on the helicity formalism have been used for all the numerical 
simulations of the channels B — > it + -n~ -V with V = K*°, K* ± ,p ± . Rigorous and detailed 
calculations of the B decay density-matrix have been carried out completely and the 
corresponding code has been already implemented in the LHCb generator code. 

• Despite the fact that the naive factorization hypothesis is very useful for weak 
hadronic B decays, this method is limited because it involves theoretical uncertainties, 
some of them being very large. 

However the physical consequences of this study are very interesting : 

• The form factor model plays important role, especially in the estimation of the 
different branching ratios which can vary by up to a factor of 2. 

• The longitudinal polarization is largely dominant, whatever the form factor model. 

• The p° — uj mixing is the main ingredient in the enhancement of direct CP violation. 

• A new way to look for CP Violation is found and it can help to develop new methods 
for measuring the angles 7 and a. 

What remains to be done is to cross-check these predictions with experimental data 
coming from the LHC experiments. 
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Figure 1: Spectrum of h ,h 00 ,h ++ . Histograms on the left correspond to the channel 

B° — > p°(u)K*° where the parameters used are: q 2 /m 2 f ) = 0.3, N^'* = 2.84, p = 
0.229, i] = 0.325 and form factors from the GH model. Histograms on the right correspond 
to the channel B + — > p°(uj)p + for the same parameters with N^- f = 2.01. 
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Figure 2: Spectrum of 38e(hij) and J^m(hij) where i ^ j. Histograms correspond to 
channel B° — > p°(u)K*° where the used parameters are: q 2 /m 2 b = 0.3, N^f = 2.84, 
p = 0.229, rj = 0.325 and form factors from the GH model. 
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Re(h _) lm(h _) 



Figure 3: Spectrum of Me(hij) and J'm{h i j) where i ^ j. Histograms correspond to 
the channel B + — > p°(u)p + where the used parameters are: q 2 /m 2 b = 0.3, N^f = 2.01, 
p = 0.229, rj = 0.325 and form factors from the GH model. 
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t 
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Figure 4: CP-violating asymmetry parameter acp{m), as a function of the tt + tt invariant 
mass in the vicinity of the oo mass region for the channel B° — > p°(uj)K*°. Parameters 
are: q 2 /m 2 b = 0.3, iV c e// = 2.84, p = 0.229,7? = 0.325. Solid triangles up and circles 
correspond to the BSW and GH form factor models respectively. 
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0.4 - 




760 765 770 775 780 785 790 795 800 



Vs (MeV) 

Figure 5: CP-violating asymmetry parameter acp{m), as a function of the ir + ir~ invariant 
mass in the vicinity of the u mass region for the channel B + — > p°(uj)p + . Parameters 
are: q 2 /m 2 b = 0.3, iV c e// = 2.01, p = 0.229, r] = 0.325. Solid triangles down and circles 
correspond to the BSW and GH form factor models respectively. 
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